Comparative analysis of specific heat of YM2B2C using nodal and two-gap models 
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The magnetic field dependence of low temperature specific heat in YNi 2 B 2 C was 
measured and analyzed using various pairing order parameters. At zero magnetic field, the 
two-gap model which has been successfully applied to MgB 2 and the point-node model, 
appear to describe the superconducting gap function of YNi 2 B 2 C better than other models 
based on the isotropic s-wave, the J-wave line nodes, or the s+g wave. The two energy gaps, 
Zlz,=2.67 meV and A 5=1.19 meV are obtained. The observed nonlinear field dependence of 
electronic specific heat coefficient, y(H)~H°' 47 , is quantitatively close to y(H)~H° 5 expected 
for nodal superconductivity or can be qualitatively explained using two-gap scenario. 
Furthermore, the positive curvature in H c2 (T) near T c is qualitatively similar to that in the 
other two-gap superconductor MgB 2 . 

PACS number(s): 74.20.Rp, 74.25.Bt, 74.25.Jb, 74.70.Dd 



Rare-earth nickel borocarbides RNi 2 B 2 C (R=Dy, Ho, 
Er, Tm, Lu, and Y) have been among the most studied 
superconductors during the past decade. The general 
interest centers around their many intriguing physical 
properties, such as the relatively high superconducting 
transition temperature T c ~\5 K, the coexistence of 
superconductivity and long-range magnetic order, and 
the reentrant superconductivity. 1 ' 2 After having 
extensive theoretical and experimental works, 3 it 
remains unclear whether they are conventional or exotic 
superconductors. 4 Particular attention has been aimed 
at determining the superconducting order parameter, 
which is thought to be essential in establishing a 
microscopic model of its superconductivity. At an early 
stage, the order parameter in YNi 2 B 2 C was considered 
to be a conventional isotropic s-wave pairing and 
mediated by conventional electron-phonon 
interactions. 5 " 8 However, recent thermal and 
spectroscopic experiments indicate the high 
anisotropy " and t/-wave pairing with line nodes. 
Furthermore, angular dependence of the thermal 
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conductivity, " specific heat, ' ' point-contact 
tunneling 17 and ultrasonic attenuation 18 all seem to 
provide evidence of point nodes in the superconducting 
gap function. Accordingly, a hybrid s + g pairing 
symmetry has been proposed for this material. 2 26 On 
the other hand, the multiband superconductivity in 
YNi 2 B 2 C has recently been proposed based on the 
observations from the temperature dependence of the 
upper critical field 27 ' 28 and the directional point-contact 
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spectroscopy. ' 

In this Letter, we present specific-heat data of 
YNi 2 B 2 C. By fitting them to various superconductivity 
models, it is found that two-gap model, which has been 
successfully applied to MgB 2 , best describes the gap 
function of superconducting YNi 2 B 2 C. 

The single crystal YNi 2 B 2 C used for the present study 
was grown by the high temperature Ni 2 B flux method. 
Details of preparation have been published elsewhere. 31 
The low temperature specific heat C(T,H) was measured 
with a He heat-pulsed thermal relaxation 



calorimeter ' in the temperature range from 0.6 to 20 
K and under magnetic fields up to 8 T. The magnetic 
field was applied in the direction perpendicular to the 
c-axis of the crystal. 
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FIG. 1. Temperature (T) dependence of specific heat (C) 
plotted as C/T vs. T 2 at various magnetic fields for YNi 2 B 2 C. 
The inset shows the deviation of the fit of 8 T data to C„(7) = 
AT^+yJ+fif+aT 5 . 

Figure 1 shows the results of calorimetric 
measurements C(T,H) of YNi 2 B 2 C as C/T vs. T 2 with H 
varying between and 8 T. Occurrence of 
superconductivity becomes evident through a specific-heat 
jump at the temperature T c for a given field. At zero field, 
the extrapolation of the data below T c to K points to a 
zero intercept, indicating a full superconducting volume of 
fraction and confirming the good sample quality. The T c 
value of 13.77 K as determined from the midpoint of 
specific-heat jump is consistent with that obtained from 
resistivity and magnetization (not shown). The 
normal-state specific heat at zero magnetic field can be 
simply described by 



CJT) = y n T+C k ,„ ic JT), (1) 

where y„T is the electronic term due to free charge 
carriers and Ci atlice (T) = pi 3 + aT 5 represents the phonon 
contribution which is assumed to be independent of the 
magnetic field. In order to achieve an optimal 
normal-state fitting to Eq. (1) for 77=8 T, the data between 
3.5 and 20 K are used but with an additional T 2 
hyperfine-contribution term, which is thought to be due to 
a very low concentration of paramagnetic centers. The 
best fitting parameters, y„ = 19.74 + 0.27 (mJ/mol K ), J3 
= 0.077 ± 0.003 (mJ/mol K 4 ) corresponding to 9 D = 533 ± 
7 K and a = 0.00018 ± 0.00001 (mJ/mol K 6 ), are fairly 
consistent with those previously reported. 7. In fact, 
these parameters are justified by the entropy balance, 

T rSC 

entropy balance S = dT at the second-order 

J j 
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superconducting-normal phase transition supporting the 
validity of the fitting above, where SC(T) = C(T,H = 0) - 
C n (T). 
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FIG. 2. A logarithmic plot of CJy„T c \s. TJT for YNi 2 B 2 C in its 
superconducting state. The solid line is the linear fit to the data 
for TJT between 4 (7=3.7 K) and 12 (7=1. IK). The left-side 
inset shows the data fitted to the BCS model. In the right-side 
inset, the almost linear relation between CJy„T and (T/T c ) 2 
indicates a nearly ^'-dependence of C e . 

The electronic specific heat in the superconducting 
state (T<T C ) is given by C e (T)= C(T)-Ci a ,,i ce (T). The 
logarithmic plot of C e (T)/y„T c vs. TJT in Fig. 2 shows a 
linear fitting of the data between TJT = 4 (7-3.7 K) and 
12 (7-1.1 K) following the relation CJT)ly n T c - exp [(—a 
TJT)] with a=0.64. The parameter a is related to the 
superconducting energy gap. Realizing that the BCS 
theory predicts C e (T)/y n T c - exp [(-1.44 TJT)] in the 
weak coupling limit, the fitted value of a is too small to 
support such a relatively high T c - 13.77 K for a single 
fully gapped BCS-like superconductor. Furthermore, data 
points in Fig. 2 deviate notably from linearity at 
temperatures higher than around 4 K. This signifies that 
YNi 2 B 2 C is not a simple BCS-like conventional 
superconductor. The poor data fitting shown in the 
left-side inset supports this argument. Significantly, 
similar observations have been observed in the two-gap 
superconductor MgB 2 . 33 ' 34 




5 10 150 5 10 15 



r(K) 

FIG. 3. Various fitting of CJT vs. 7 using (a) line-node, (b) 
point-node, (c) s+g wave, and (d) two-gap models. The 
quality of each fitting is shown in the respective inset, where DF 
represents the difference between the calculated values and 
actual data. 

To further explore the possibility of other 
unconventional superconducting order parameters in 
YNi 2 B 2 C, the data are fitted to line-node, point-node, s+g 
wave and two-gap models, as shown in Figs. 3(a), (b), (c) 
and (d), respectively. The equation used in evaluating the 
superconducting-state electronic specific heat C e for these 
models is 

C e =2N(0)J3k~ jdijddsmff J -^-(E 1 +- fi—)de' ^ 
^ n o o -nm a 2 dp 

where N(0) is the density of states at the Fermi surface, 
f}=\lkT, £=(£ 2 +A 2 )" 2 , f={\+e^) \ A=A the 
superconducting energy gap for isotropic s wave, 
A=A o cos«0 for line nodes (to simplify, n=2 for rf-wave), 
A=A sinn<? for point nodes and A=A /2(l-sin 4 6' cos4$) 
for s+g wave. For simplicity, the Fermi surface is 
assumed to be spherical except that, for the case of the 
line nodes, a two dimensional Fermi surface is applied. 
The deviation of the fit (DF) in the inset indicates clearly 
that the line-node model can not satisfactorily describe the 
data as shown in the Fig. 3(a). The fit of data to the 
point-node model is quite good at low temperatures, 
unless it becomes a little ambiguous at the high 
temperature region as shown in the Fig. 3(b). This can be 
related to the fact that C e in Eq. (2) is proportional to T 3 
only at low temperatures in the point-node case, while the 
observed Independence of C e prevails almost throughout 



the entire superconducting state (T < T c ) as shown in the 
inset of Fig. 2 and also in several previous reports. 9 ' ' 3 
However, considering the experimental resolution, the 
point-node symmetry cannot be totally ruled out from the 
present specific-heat study. On the other hand, Eq. (2) 
gives C e ~ I 2 at low temperatures for the s+g wave case 
(as addressed in Ref. 26), thus the unacceptable fit is seen 
in the Fig. 3(c). Furthermore, from the analysis of the 
H=8 T data, the concentration of paramagnetic impurity is 
less than 10 3 . Therefore, even considering the scenario of 
impurity scattering, 26 the observed T -dependence is 
difficult to reconcile with the s+g wave model. Other than 
these three possibilities, the two-gap model, which has 
been successfully applied to MgB 2 superconductor, 3 ' 3 
was employed to further analyze the present data. In this 
two-gap scenario, two distinct gaps 2 A i/k B T c =4.5 and 2 
A <A fi 7>2.0 with 71% and 29% in relative weight, 
respectively, are introduced. The excellent fit to the data is 
shown in Fig. 3(d). Consequently, the corresponding large 
gap A i=l. 67 meV and the small gap Zl s =1.19 meV can 
be obtained. It is noted that the derived gap values are 
fairly consistent with the recent result from point contact 
spectroscopy, 30 where an anisotropic (or a multiband band) 
gap withZl max =2A meV and A „„■„ =1.5 meV was revealed. 
Conceptually, the point-node and two-gap models are 
incompatible. However, for a highly anisotropic case it is 
not easy to experimentally distinguish between them. 
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FIG. 4. Magnetic field dependence of electronic specific heat 
coefficient y(H) and the solid line representing y(H) ~ H° is 
the best fit. The inset shows the determination of y(H) from the 
linear extrapolation of data for each field below 2 K in Fig. 1. 

An alternative approach to study the superconducting order 
parameter is through the vortex excitations in the mixed state. 
Figure 4 shows the linear coefficient of electronic specific heat 
y(H) as a function of the applied magnetic fields. Each point was 
obtained from the linear extrapolation of data below 2 K for a 
given field as shown in the inset. The best fit yields y(H) ~ H 047 . 
Apparently, y(H) follows an H" dependence with a being very 
close to 0.5 as reported previously. 916 ' 20 - 21 Such a y(H) ~ H lh 
relation was used to support that YNi 2 B 2 C is possibly a nodal 
superconductor 20 " 22 similar to a nodal rf-wave cuprate 
superconductor 32 ' 36 as predicted by Volovik. 37 However, the 
nonlinear H dependence of y{H) has also been observed and 
successfully argued to be an intrinsic property of the two-gap 
superconductor MgB 2 . 33 ' 34 ' 38,39 
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FIG. 5. Temperature dependence of upper critical field H c2 (T) 
of YNi 2 B 2 C obtained from Fig. 1. The crosses representing data 
from magnetization measurements provide a consistency check. 
The dashed line simply connects the data points, while the solid 
line is calculated from the WHH theory. Inset shows the 
determination of error bars. 

Figure 5 shows the temperature dependence of the upper 
critical field H c2 (T) of YNi 2 B 2 C obtained from Fig. 1, 
where the solid line is based on the 
Werthamer-Helfand-Hohenberg (WHH) theory. The 
salient features are the positive curvature near T c and all H c2 
values being smaller than the theoretical prediction. 
Meanwhile, the H c2 (T) behavior is consistent with that 
measured resistively on YNi 2 B 2 C and LuNi 2 B 2 C, 28 from 
which Shulga et al. have successfully calculated such an 
H c2 (T) behavior based on the two-band model of 
superconductivity. 28 Just as convincing, similar 
observations have also been proposed to reflect the intrinsic 
property of two-gap superconductivity in MgB 2 . 40 

In summary, the magnetic-field dependence of low 
temperature specific heat of YNi 2 B 2 C has been analyzed 
using various superconducting symmetry order parameters. 
Several critical findings include: (1) At zero magnetic field, 
the superconducting gap function of YNi 2 B 2 C is better 
described by the point-node and two-gap models than other 
models based on single isotropic s wave, line nodes or s+g 
wave. (2) The two gap values, Zl/,=2.67 meV and 
Zl 5= 1 . 1 9 meV obtained from the data fitting to the two-gap 
model are consistent with recent reported multiband gaps 
determined from point-contact spectroscopy. The observed 
nonlinear relation y{H) ~ H ' 47 in the mixed state can be 
explained by either nodal or two-gap superconductivity. (3) 
The positive curvature observed in H c2 (T) near T c is similar 
to that in MgB 2 . These findings provide calorimetric 
evidence to support that YNi 2 B 2 C, might be a point-node or 
likely as MgB 2 a two-gap superconductor. 
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